Relativistic corrections to isotope shift in light ions. 
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We calculate isotope mass shift for several light ions using Dirac wave functions and mass shift 
operator with relativistic corrections of the order of (aZ) 2 . Calculated relativistic corrections to the 
specific mass shift vary from a fraction of a percent for Carbon, to 2% for Magnesium. Relativistic 
corrections to the normal mass shift are typically smaller. Interestingly, the final relativistic mass 
shifts for the levels of one multiplet appear to be even closer than for non-relativistic operator. 
That can be important for the astrophysical search for possible a-variation, where isotope shift is 
a source of important systematic error. Our calculations show that for levels of the same multiplet 
this systematics is negligible and they can be used as probes for a-variation. 

PACS numbers: 06.20.Jr, 31.30.Gs, 31.25. Jf 



I. INTRODUCTION 

Modern theories (such as the string theory and M- 
theory) predict temporal and spatial variations of funda- 
mental physical constants. Some recent studies of quasar 
absorbtion spectra indicate that the fine structure con- 
stant a = e 2 /he, where e is the charge of the electron 
and h and c are the reduced Planck constant and the 
speed of light respectively, could have changed during 
the evolution of the Universe (see [l[ and the references 
therein). This result was not confirmed by other groups 
@; HI, EL IE @j an d hence new experiments are required. 
Laboratory experiments are rapidly increasing their sen- 
sitivity to a-variation (see and references therein), 
but are currently slightly less sensitive than astrophysi- 
cal ones, where the time scale is 10 10 times longer. Fi- 
nally, there may be some evidence for varying a from the 
natural nuclear reactor in Oklo, which operated about 2 
billion years ago [§]. 

Astrophysical studies are based on the fact that 
atomic transition frequencies depend on the parameter 
x = (a/ao) 2 — 1, where ao ~ 1/137 is the laboratory 
value of a. Thus, one can look for a-variation by com- 
paring the frequencies of atomic lines in the spectra of 
quasars with their laboratory values, which correspond 
to x — 0. One of the major systematic effects, which can 
imitate a- variation, is the isotope shift (IS). It had been 
shown previously [l(| HH, [H, El that typical IS is of the 
same order of magnitude as the effect observed in [l| . 

In [14| . a method to separate the effect of varying a 
from the IS was proposed. This method is based on build- 
ing special combinations of atomic frequencies, some of 
which are insensitive to both effects ("anchors"), while 
the other are independent on the IS, but sensitive to a- 
variations ( "probes" ) . One could then use the anchors to 
retrieve the actual value of the redshift of the spectrum, 
and then use the probes to evaluate the size of the varia- 
tion Aa/a. In addition, one can form the probes that are 



'Electronic address: vakorol@mail.ru 



sensitive to the IS, but not to the variation of a. These 
probes could provide information about isotopic abun- 
dances in the early Universe. Such information should 
also be of interest to the astrophysicists, because isotopic 
evolution is tightly linked with the general evolution of 
the Universe. 

Described method of anchors and probes requires 
precise calculations of atomic relativistic coefficients q, 
which determine the dependence of transition frequen- 
cies on x, and the IS coefficients fcis- Factors q strongly 
depend on the total angular momentum J, which distin- 
guishes levels of the same multiplet. Relativistic correc- 
tions to the coefficients k\§ are known to be small and 
their dependence on J is much weaker. Nevertheless in 
order to separate two effects reliably, one needs to accu- 
rately account for relativistic corrections to IS. For light 
atoms IS is dominated by the mass shift (MS). Usually 
MS calculations are done with Dirac wave functions, but 
using non-relativistic form of MS operator. In this pa- 
per we report results of the MS calculations for 12 light 
atoms and ions. We use the MS operator, which includes 
relativistic corrections of the order of {otZ) 2 fl5l. [lfjl. \vf\ . 
It is shown that relativistic corrections to the operator 
partly compensate corrections to the wave functions and 
final relativistic corrections are smaller than in previous 
calculations. That may help to distinguish between a- 
variation and isotope effects in astrophysics. 



II. METHOD 
A. Theory of isotope shift 

Total IS consists of MS and field shift (FS). The former 
accounts for the finite mass of the nucleus, and the latter 
accounts for the finite size of the nucleus. In the non- 
relativistic theory the MS is further divided into normal 
mass shift (NMS), and specific mass shift (SMS). The 
overall IS in the frequency of any atomic transition of 
the isotope with mass number A' with respect to the 
isotope with mass number A can be written as (see, for 
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example 



' A = (fcNMS + fcsMs) 



1 X +FS(r 2 f' A , (1) 
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where (r 2 ) is the nuclear mean square radius, F is the FS 
coefficient and &ms = ^nms + ^SMS is the MS coefficient. 
In the non-relativistic theory NMS is described by one- 
electron operator, 



(nr) 
NMS 



1 

2M 



2^ Pi 



(2) 



where pi is momentum of the i-th electron. Because of 
the virial theorem Anms has the simple form: 



k 



(nr) 
NMS 



1823' 



(3) 



where the number 1823 refers to the ratio of the atomic 
mass unit (amu) to the electron mass. SMS is described 
by a two-electron operator: 



H, 



(nr) 
SMS 



Pk 



(4) 



in- 



consistent relativistic theory of MS can be formulated 
only within quantum electrodynamics in a form of ex- 
pansion in aZ. It was shown in Refs. [HI, [l7| that 
the first two terms of the expansion in aZ lead to the 
following relativistic MS Hamiltonian: 



(r) 
MS 



(5) 



where is the Dirac matrix of the i-th electron and 
f = r jr. As in the non-relativistic theory, one can divide 
Hamiltonian J5|) into the one- and two-electron terms: 
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Non-relativistic expression ([3]) is not applicable to oper- 
ator ([6]) and now /cnms has to be calculated on a same 
footing as Asms- 

As pointed out above, for the light atoms and ions FS 
is much smaller than MS. Still for the comparison with 
high precision experiments one may need to account for 
FS. Below we focus on calculating MS, but in the final 
Table IVl we present FS factors for some of the transitions. 



B. Electron correlations 

IS is very sensitive to electron correlations. High accu- 
racy calculations must account not only for correlations 



between valence electrons, but also for core-valence cor- 
relations. Here we treat both types of correlations within 
CI+MBPT method [H]. In this method the configura- 
tion interaction (CI) calculation is done for valence elec- 
trons using effective Hamiltonian H c g, which is formed 
within the second order many body perturbation theory 
(MBPT) in the residual core-valence interaction. All our 
calculations are done with modified Dirac- Fock code [2(| , 
CI code (U, and MBPT code 0. 

CI+MBPT method is easily reformulated for the IS 
calculations within the finite-field approximation. In this 
approximation the IS operator His is added to the many- 
particle Hamiltonian H with an arbitrary coefficient A: 



H + XH 



is- 



(8) 



The eigenvalue problem for Hamiltonian |(5J| is solved for 
+A and for —A. Then, the IS correction to the energy is 
recovered as: 



AE ls = 



E+x - E. 
2A 



(9) 



All calculations are done independently for the field, 
normal, and specific parts of the IS. Parameter A is cho- 
sen from the considerations of the numerical stability and 
smallness of the the non-linear terms. In the CI+MBPT 
calculations Hamiltonian © is used to construct H\ e g. 
That means that H\ is used on all stages of calculation 
starting from solving Dirac-Fock equations for atomic 
core. 

CI+MBPT method provides the solution of the eigen- 
value problem for valence electrons H c s^ = E^>, which 
is an approximation to all-electron problem H^> = E^S . 
The latter equation gives an exact solution for the given 
basis set. Unfortunately this equation can be solved only 
for atoms with few electrons and only for very short ba- 
sis sets. Nevertheless, on a short basis set we can com- 
pare all-electron and CI+MBPT calculations with some 
simple calculation (for example, with one-configurational 
approximation) and designate corresponding corrections 
as A short and A^f rt . Then we repeat CI+MBPT calcu- 
lation on a much longer basis set and find new correction 
Af ff ng . Although we can not make all-electron calculation 
for the long basis set, we can use the short basis set cal- 
culation to shift the central point of the final CI+MBPT 
calculation and estimate its error: 



<$shift 
^crr 



short 
^shift 



Ashort 
eff 



A lo " s 

cti 
^short 
err 



(io) 
(ii) 



We will use Eq. (|I0|) where possible to improve 
CI+MBPT results and Eq. (TT]) to estimate theoretical 
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III. RESULTS OF THE CALCULATIONS 
A. Lithium-like ions 

Lithium is the most simple atom in our consideration, 
since it has only one core orbital Is and one valence elec- 
tron. The simplicity of this system allows one to perform 
all-electron (full) CI calculation for all three electrons. 
Such calculation provides an exact solution on a given 
one-electron basis set and allows to improve our results 
and to control their accuracy using Eqs. (fTOlfTTj) . The ba- 
sis set consists of the orbitals: 1-9 s, 2-9p, 3-8d, and 4-8/ 
(basis set [9sp8df]). CI+MBPT calculation is then done 
with two longer basis sets: valence CI space is defined 
by the basis set [13spl2d/] and all intermediate MBPT 
summations run over basis set [17spl8c?19/20<721ft,] (see 
for details). 

In Table Q] we compare our MS calculations with the 
experiment. For Li the relative size of FS is about 1CP 5 
and we neglect it. 



TABLE I: Isotope shifts (in MHz) in the spectrum of 6 Li, 
with respect to 7 Li. In the second column result of our 
CI+MBPT calculation is given. Third column presents cor- 
rected values using Eq. (|10[) and the errors are estimated with 
the help of Eq. JTT). 

Transition Theory Exper. 

CI+MBPT Final 



2Si/ 2 


-> 2pi /2 


10346 


10608(300) 


10534.26(13)" 
10532.9(6)'' 


















10534.3(3) c 










10533.13(15)'' 
10533.3(5) 6 
10539.9(1.2) c 
10534.93(15) d 


2Sl/2 


-» 2p 3 /2 


10346 


10607(300) 


2Sl/2 


-> 3pi /2 


14162 


14354(200) 


14470(450) e 


2Si/ 2 


-> 3p 3/2 


14162 


14398(250) 


14470(450) c 


2si/ 2 


-> 3d 3/2 


13194 


13450(300) 


13314(6/ 
13312(4) 9 


2si/ 2 


— * 3d 5 / 2 


13194 


13450(300) 


13314(6/ 
13312(4) s 


2si/ 2 


-* 4s 1/2 


14560 


14750(200) 


14656(e) 7 
14661(14) 9 



a Walls et al.\%& 
'Sansonetti et q/.[24| 
c Windholz and Umfer [2I 
d Scherf et aZ.[HJ 
e Mariella [13 
^Lorenzen and Niemax [28 
SKowalski et aZ.fH 



We see that most experimental results for IS in Li are 
much more accurate than our calculations. Corrections 
(|10[) lead to significant improvement of the agreement be- 
tween calculations and experiment. The error estimate 
(fTTj) appears to be quite reliable here as most calculated 
CI+MBPT values and all corrected values agree with ex- 
perimental data within the estimated errors. We con- 
clude that Eqs. (fTUI [TTj) work nicely for Lithium. On the 



other hand, even for such a simple system as Li, our the- 
oretical error for MS appears to be about 2%, while the 
typical error for the frequencies is only about 0.1%. That 
shows the extreme sensitivity of the MS to correlations 
and significantly more important role of the high order 
MBPT corrections, which are neglected here. 

Calculations of other Li-like ions are done in a similar 
way using basis sets of the same length. Corresponding 
results for C iv, N v and O vi are given in the final Ta- 
ble [V] There we also use Eq. (jlip to estimate theoretical 
accuracy. 



B. B II, C III, and C II 

We treat Be-like ions B 11 and C in as two-electron ions 
with Is 2 core and the ground state ^So [2s 2 ] . Since the 
full 4-electron CI calculation is more complex than in the 
case of 3-electron Li, a shorter basis set [5spAdf] is used 
here for such calculation. Valence 2-electron CI calcu- 
lations were done with [10sp9df] and MBPT corrections 
were calculated using [16spl7dl8/19<?20/i] basis set. 

There are many experimental results for B II. Table HT1 
presents the comparison of calculated and experimental 
IS for B 11. Unlike the case of lithium, here for most 
transitions the estimated theoretical error is smaller than 
experimental one. We see that generally there is good 
agreement between theory and experiment. Our calcula- 
tion of IS in transition ]P? — > 1 Z?2 is in good agreement 
with the measurement [30] ? but is 3<r away from the re- 
sult of Ref. [H|. Note that for the transition l P£ -> % 
the MBPT correction for the long basis set appears to 
be anomalously large. That leads to the much larger 
theoretical error, than for other transitions. 

C 11 is the only B-like ion which is interesting for astro- 
physics. It is also of a particular interest to us because 
the experiment [33j hints at approximately 2% difference 
between the shifts for the relativistic doublet. If this is 
confirmed, it will mean that relativistic corrections to IS 
are quite large even for light ions and should be included 
in accurate calculations. That will also make isotope ef- 
fects more dangerous for the astrophysical searches for 
a- variation. 

The ground state of C 11 is 2 P°/ 2 [2s 2 2p!/ 2 ] . It has five 
electrons and full CI becomes impractical even for the 
short basis set. That makes it more difficult to deter- 
mine theoretical error. Using our calculations for Li-like 
and Be-like ions we estimate theoretical error for C 11 to 
be about 1 - 2%. On the other hand, with relativistic MS 
operators ^ and the relative size of the shifts within 
one multiplet should be also accurate to few percent. Our 
results shown in Table UlTl indicate that the difference be- 
tween the shifts in this doublet must be roughly two or- 
ders of magnitude smaller than hinted by the experiment 
[33l ]. In fact, as we discuss below in more detail, using rel- 
ativistic operator ([5]) brings MS shifts for the levels of one 
multiplet closer to each other than with non-relativistic 
operators (J2J) and 
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TABLE II: Isotope shifts (in MHz) for n B with respect to 
10 B in the spectrum of B II. 



Transition 






IS 






A (A) 


Theory 


Exper. 


'Softs 2 ] -> L P?[2s2p] 


1362 


21303(66) 


21000(3000)" 
21000(3000)' 




^ 3 P 2 [2p 2 ] 






3 P?[2s2p] 


1624 


22367(170) 


21600(3600)" 
24000(3600)' 




-> 3 Pi[2p 2 ] 






3 PS[2s2p] 


1624 


22369(170) 


21600(3600)" 
21600(3600)' 




-> 3 Pi[2p 2 ] 






3 P?[2s2p] 


1624 


22370(170) 


21600(6900)" 
21600(6900)' 










3 Pi[2s2p] 


^ 2 P 2 [2p 2 ] 


1624 


22369(170) 


21600(4500)" 
18300(4500)' 




-» 3 P)[2p 2 ] 






3 P?[2s2p] 


1624 


22371(170) 


21600(3600)" 
21600(3600)' 




^ 3 Pi[2p 2 ] 






3 P 2 °[2s2p] 


1624 


22372(170) 


21600(3600)" 
21600(3600)' 
12600(1200)' 
26600(100) ' 
26300(120) c 










1 A°[2s2p] 


1 5o[2p 2 ] 


1843 


12689(580) 


'PiP&Pl 


3451 


26663(55) 



"Jonsson et al.\32 
'Litzen et al. [3y] 
c Vinti [H 



TABLE IV: Comparison with experiment for IS in Na I and 
Mg I, II. The FS is neglected (see text). All numbers are in 
MHz. 



Transition 


A, A 


MS 


IS (expt.) 




^ 3 Na I 






3s —> 3pi /2 


5896 


775.8 


758.5(7)" 
756.9(1.9)' 
757.72(24) c 


3s -> 3p 3/2 


5890 


776.5 


24- 


*°Mg ii 






3s -> 3p 3/2 

24- 


2796 


3086.3 


3050(100)^ 



^opa 2 ] - 
'Sops 2 ] - 
3 Po°[3s3p] 
3 P?[3s3 P ] 
3 P 2 °[3s3p] 
3 Pi°[3s3p] 
3 P 2 °[3s3p] 
^ops 2 ] - 



3 Pi°[3s3p] 4571 
'PTpsSp] 2852 

^ 3 Si[3s4s] 5167 

^ 3 Si[3s4s] 5173 

-> 3 Si[3s4s] 

+ 3 D 2 [3s3d] 

+ 3 P 3 [3s3d] 

'pr&sM 

2T 



5184 
3832 
3838 
2026 



2686.9 
1425.3 
-367.9 
-369.9 
-373.1 
58.7 
55.8 
2975.0 



^ops 2 ] - 
%ps 2 ] - 
3 Pfps3p] 
^ops 2 ] - 



'Mg i 
3 P°ps3p] 4571 
'PipsSp] 2852 
^ 3 5ips4s] 5184 
1 P 1 °ps4p] 2026 



1397.2 
741.2 
-194.0 
1547.0 



2683.2(0) e 
1415.3(5.0) / 
-395.7(6.0) s 
-391.2(4.5) 9 
-392.7(7.5) 9 
60.6(3.0)'" 
58.2(3.6) h 
2918. 9(1. 2Y 



1405. 2(0. l) e 
743.8(3.0/ 
-203.9(7.5) 9 
1526.4(1.2)' 



TABLE III: IS for 13 C in respect to 12 C in the doublet line 
2 S 1/2 [2s2p 2 ] -»• 2 P^[2s 2 3p] of C II. Here we use relativistic MS 
operators ([§]) and ©• Calculation of Berengut et al. [34|] was 
done with non-relativistic operators. 



Transition 


A, A 


IS (MHz) 
This work Theory [34] 


Exper. [33] 


S>l/2 ^3/2 
1q 2 ryo 
0\/2 > ^1/2 


2836.707 
2837.605 


-18185.0 -18500 
-18185.9 -18500 


-18350(60) 
-18680(90) 



C. Sodium-like ions and Mg I 

There are several Na-like ions which are important for 
astrophysics, namely Na I, Mg II, Al III, and Si iv. For 
Na I and Mg II there are experimental data available 
and we compare calculated MS with the experiment in 
Table IIVI Using nuclear radii from Ref. [35| we have 
estimated FS for 3s — > 3p transitions in 22_23 Na I and 
24 ~ 26 Mg II to be about -2 and +16 MHz, respectively. 
This is significantly smaller than the uncertainty in our 
MS calculation. 

In these calculations for Na I, Al in and Si iv 
we use CI basis set [14spl3c?f] and MBPT basis set 
[19spl8dl9/20g21/i]. Our IS results for Si iv are within 
1% agreement with calculation of Berengut et al. [3t| . 
while the differences for Mg I are larger and can reach 
several percent. 

Neutral magnesium is one of the most well studied 
two-electron atoms. The data from several experimen- 
tal works on IS in Mg i,ii are presented in Table HVl and 
compared with theoretical IS calculated on the basis sets 
[12spd9/] and [17spdf9gh]. 



"Pescht et al.pi^ 
'Huber et al.)M 
c Gangrsky et al. [3|1 
d Drullinger et a^.l40l 
e Sterr et al.]£j& 
-fBeverini et a/.li^l 
9 Hallstadius & Hansen 
h Hallstadius [4^1 
'Hanneman et al. 45] 



IV. DISCUSSION 

Our final results for MS in astrophysically important 
transitions are listed in Table [V] Where available we 
also give FS parameters and results of other calculations 
of SMS and relativistic g-factors. The latter are known 
to be significantly different for the lines of one multiplet 
because the splittings within the multiplet are caused 
by relativistic spin-orbit interaction and scales as (aZ) 2 . 
That makes these splittings natural candidates as probes 
for cv-variation. However, J-dependence of IS may in- 
troduce some systematic errors as isotope abundances in 
the Universe may vary significantly. Here we use rela- 
tivistic MS operators and show that J-dependence for 
MS is smaller than one might expect from the simple 
considerations and from experiment [33J. That means 
that multiplet splittings are indeed good probes for cv- 
variation. 

In several cases we also made SMS calculations using 
non-relativistic operator (|4]) (see Table IV)) . Comparison 
with our relativistic values indicate that relativistic cor- 
rections grow from a fraction of a percent for Carbon ions 
to about 1 - 2% for Mg i,ii. It is interesting though, that 
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TABLE V: Isotope shift parameters for astrophysically important transitions in Li-, Be-, B-, and Na-like ions, and in Mg I. The 
experimental values for the transition wavelengths were taken from the NIST atomic spectra database (http://physics.nist.gov). 
The SMS and NMS parameters are given in GHz-amu, the FS parameter F is given in MHz/fm 2 . Non-relativistic values for 
/csms and /csms were obtained with operator Q and scaling @ respectively. 



Transition 


A, A 


non-rel. 


fcsMS 

rel. other works 


non-rel. 


/cnms 

rel. 


F 


q 










C IV 










2S1/2 


->■ 2p 1/2 


1551 


-4509 


-4498(45) -4511(23)" 


-1060.5 


-1061.1 


-174 


104 d 


2S1/2 


-»■ 2p 3/2 


1548 


-4502 


-4499(45) -4504(23)" 


-1062.3 


-1061.2 


-173 


232 d 


'Sops 2 ] 

'Sops 2 } 
'Sops 2 ] 


-» A PS [2s2p] 
-» 3 P! [2 S 2p] 
-» 3 P 2 °[2s2p] 


1910 
1909 
1907 
977 


-3461 
-3459 
-3456 
-2778 


C m 

-3453(7) -3473" 
-3454(7) -3472" 
-3454(7) -3468" 
-2774(80) -2790" 


-861.2 
-861.6 
-862.6 
-1683.3 


-860.3 
-860.2 
-860.1 
-1677.3 




74" 
108" 
178" 
165" 













C ii 










2 P° /2 [2 S 


2 2p] 




2 S 1/2 [2s2p 2 } 


1036 


-1295 


-1321" 


-1586.9 


-1591.1 


168" 


2 P° /2 [2s 


2 2p] 




2 D 3/2 [2s2p 2 ] 


1335 


-2636 


-2671" 


-1232.3 


-1230.3 


178" 


2 P° /2 [2 S 


2 2p] 




2 D 5/2 [2s2p 2 ] 


1334 


-2637 


-2672" 


-1232.3 


-1230.4 


181" 


2 P° /2 [2 S 


2 2p] 




4 P 1/2 [2s2p 2 ] 
4 P 3/2 [2 S 2p 2 ] 


2325 


-2971 


-2960" 


-707.2 


-700.8 


132" 


2 P° /2 [2 S 


2 2p] 




2324 


-2972 


-2958" 


-707.6 


-700.8 


158" 


2 P° /2 [2 S 


2 2p] 




4 P 5/2 [2s2p 2 ] 


2322 


-2973 


-2956" 


-708.1 


-700.6 


202" 











N v 












2si/ 2 


-> 2p 1/2 


1243 


-7123 


-7099(70) 




-1323.3 


-1324.7 


-417 


196 e 


2si/ 2 


-> 2p 3/2 


1239 


-7109 


-7103(70) 




-1327.6 


-1325.7 


-386 


488 e 










O vi 












2S1/2 


-» 2p 1/2 


1038 


-10312 


-10266(100) 




-1585.0 


-1587.0 


-719 


340 d 


2S1/2 


-»■ 2p 3/2 


1032 


-10283 


-10273(100) 




-1593.7 


-1590.5 


-719 


872 d 










Nai 












3Si/2 


-* 3pi /2 


5896 




-116 


-108(24) 6 
-97 c 


-278.9 


-278.3 


-33 


45 d 


3Sl/2 


-> 3p 3/2 


5890 




-116 


-107(24) b 


-279.2 


-278.5 


-33 


63 d 








-97 c 


















Mg ii 












3*1/2 


-» 3p 1/2 


2803 


-387 


-378 


-373(12) 6 


-586.6 


-586.7 


-127 


120 c 








-362 c 










3Sl/2 


-> 3p 3/2 


2796 


-381 


-378 


-373(6)'' 
-361 c 


-588.1 


-587.3 


-127 


211 c 










Mg i 












'Softs 2 ] 
'Softs 2 ) 


-> 'P?fts4p] 


2026 


-102 


-101 


-108 / 


-811.6 


-816.9 


-175 


87 f 


-» 'P?fts3p] 


2852 


129 


130 


134 ; 


-576.5 


-575.3 


-175 


94 / 










Al in 












3*1/2 


-> 3pi /2 


1863 




-835 




-882.9 


-882.6 


-254 


216 c 


3S1/2 


-»■ 3p 3/2 


1855 




-834 




-886.7 


-884.6 


-254 


464 c 










Si iv 












3Si/ 2 


-»■ 3pi /2 


1403 




-1510 - 


1535(11)' 
-1505(7) 6 


-1172.4 


-1171.1 


-497 


362 e 


3Sl/2 


-»■ 3p 3/2 


1394 




-1507 


-1180.0 


-1175.6 


-480 


766 e 



"Berengut et al. 34| 
6 Berengut et al. 3a | 
c Safronova & Johnson 
d Berengut et al. 4al 
e Dzuba et aZ.ggl 
^Berengut et al. 5(|] 
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relativistic MS values for one multiplet tend to be closer 
to each other than non-relativistic values. The difference 
between final MS shifts within one multiplet appear to be 
much smaller than relativistic corrections themselves. It 
should be mentioned that non-relativistic SMS presented 
in Table [V] was obtained by applying operator |4]) only 
to the upper components of the Dirac functions. If this 
operator is applied to both components, the result tends 
to be noticeably farther from the correct relativistic value 

It should be mentioned that relativistic calculations 
of NMS require special care. Indeed, even in the non- 
relativistic limit calculated NMS will correspond to the 
scaling ([3]) for theoretical, rather then experimental fre- 
quency. For light many-electron atoms the difference 
between experimental and calculated frequencies can be 
larger than actual relativistic effects. To account for that 
one can multiply calculated relativistic NMS by the ratio 
^exper/wtheor- The values presented in Table [V] include 
this correction. 

Using q- factors from Table fVl we can also form number 
of probes for isotope abundances, which are insensitive to 
a- variations. There are several such probes for Carbon 
ions. The following probe for C IV: 



These probes have &MS = 1810 GHz-amu, 776 GHz-amu 
and 1593 GHz-amu, respectively. Finally, one can form 
two probes for Nitrogen and Oxygen: 



pis 



pis 



0.71w(2a 1/2 - 2 Pl/2 ) Nv (16) 
0.29w(2s 1/2 -> 2p 3/2 ) Nv = 33719.52 cm- 1 

(&ms = 3537 GHz-amu), 
0.72^(2 Sl/2 ->2p 1/2 ) Ovi (17) 

0.28w(2s 1/2 -> 2p 3/2 ) Qvi = 42255.9 cm" 1 
(&ms = 5213 GHz-amu). 



V. CONCLUSIONS 



Pl s = 0.69w(2si /2 
- 0.31cj(2si/ 2 



'Civ 



2pi/2j 

2 P3/ 2 ) CIV 



(12) 



24470.71 cm" 1 , 



has the MS parameter /cms = 2112 GHz-amu. More 
probes can be constructed from C in lines: 



pIS 


= 0.71w( 1 S'o - 


-> 3 < 


'cm 




(13) 




- 0.29w( 1 S* - 


- %° ; 


'cm 


= 21970.95 cm" 1 , 




pIS 


= 0.59a; (^o - 




'cm 




(14) 




- 0.41w( 1 S'o - 




'cm 


= 9416.36 cm" 1 , 




pIS 

■M 


= OMlu^So - 




'cm 




(15) 




- 0.31w( 1 S* - 




'cm 


= 4404.14 cm -1 , 





We have used the method of the effective Hamiltonian 
to calculate IS for astrophysically important transitions 
in ions with several valence electrons. We used relativis- 
tic MS operators and found out that relativistic correc- 
tions for light ions are not suppressed numerically and are 
of the order of (aZ) 2 , i.e. 1% for Z « 10. On the other 
hand, the difference between MS values for one multiplet 
appears to be at least one order of magnitude smaller. 
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